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1 .0 SUMMARY 

This document provides a set of preliminary entry guidance and autopilot software 
formulations for use in the Mission Control Center (MCC) entry processor. These 
software formulations meet all level B requirements, as specified in reference 1. 


2.0 INTRODUCTION 

This internal note presents the level C software formulation requirements for 
the entry guidance and the simplified autopilot that will be used by the MCC 
entry processor. The entry guidance logic is based on the Orbiter avionics 
entry guidance software, as described in reference 2. This MCC requirements doc- 
ument contains a definition of coordinate sy stems (sec. 3.2), a list of parame- 
ter definitions for the software fo rmulations (sec. 3 - 3 ) , a desc ription of the 
entry ~~guTcTance deta IT ed fonnul ation requirements (sec. 3 .4) , a description of 
t^he detailed autop ilot for mulation requirements (sec. 3 -5 ) , a description of the 
targe ting r outine (sec. 3.6), and a set of formulation flow charts (appendixes 
X~through C ) . 


3.0 SOFTWARE FORMULATION REQUIREMENTS 


3.1 INTRODUCTION 

The entry guidance system is the source of the bank-angle and angle-of-attack 
commands used to control the entry trajectory. The entry guidance can be called 
by one of two means: either in the normal mode to generate a complete entry 

trajectory, or as a part of the iterative targeting mode in the entry target gen- 
eration (ETG) subphase of the deorbit processor. This logic will be discussed 
in section 3.4. 

The autopilot generates the Orbiter attitude response to the entry guidance 
commands. This is accomplished by means of a simple phase plane in the bank- 
angle and angle-of-attack axis. The sideslip angle (&) is always assumed to be 
zero. This logic is discussed in section 3.5. 

The targeting logic generates the range and heading information to the targeted 
runway and is used by the entry guidance and the deorbit processor. This logic 
is described in section 3-6. 


3.2 COORDINATE SYSTEMS 

Two basic coordinate systems are assumed by the software formulations described 
in this document. The state vector is assumed by the targeting routine (EGRT) 
to be in the Greenwich true-of-date system, as defined in figure 3.2-1. The 
runway coordinate system is defined in figure 3-2-2 and is used by EGRT. The 
bank-angle and angle-of-attack commands generated by the entry guidance and 
executed by the autopilot are attitudes defined with respect to the Earth-rela- 
tive velocity vector. 
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Name: Greenwich true of date (geographic). 

Origin: The center of the Earth. 

Orientation: The Xg-Yg plane is the Earth's true-of-date equator. 

The Zg-axis is directed along the Earth's true-of-date 

rotational axis and is positive north. 

The +Xg-axis is directed toward the prime meridian. 

The Yg-axis completes a right-handed system. 

Characteristics: Rotating, right-handed, Cartesian. Velocity vectors expressed 

in this system are relative to a rotating reference frame 
fixed to the Earth, whose rotation rates are expressed 
relative to the Aries-mean-of-1950 system. 


Figure 3-2-1.- Greenwich true of date (geographic). 
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Name: Runway coordinate system. 

Origin: Runway center at approach threshold. 

Characteristics: Rotating, Earth referenced. 

Description: Z R -axis is normal to the ellipsoid model through the runway 

centerline at the approach threshold and positive toward the 
center of the Earth. X R -axis is perpendicular to the Z R -axis 
and lies in a plane containing the Z R -axis and the runway 
centerline (positive in the direction of landing). 

Y R -axis completes the right-handed system. 

hi is the runway azimuth, measured in the X R -Y R plane from 
true north to the +X R -axis (positive clockwise). 


I Figure 3.2-2.- Runway coordinates. 

'*! 
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3.3 PARAMETER LISTS AND DEFINITIONS 

A complete list of all parameters used in the entry guidance, autopilot, and 
targeting routine is presented with appropriate definitions. Tables 3-3-1 and 

3.3- 2 present the input and output data for the entry guidance. Table 3.3-3 
presents a list of internal parameters for entry guidance. Tables 3-3-4, 

3.3- 5, and 3.3-6 present the same data for the autopilot, and tables 3-3-7, 

3.3- 8, and 3.3-9 present the data for the targeting routine. The data in these 
tables represent STS-1 cycle 3 trajectory data. This document will not be 
updated for constant changes. These changes are documented in system parameter 
lists by the Ground Data Systems Division (3DSD). 


TABLE 3.3-1.- ENTRY GUIDANCE INPUT DATA 
<a) Input parameters 
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Symbol 

Description 

Source 

Unit 

ALPHA 

Angle of attack 

Autopilot 

deg 

DELAZ 

Current azimuth error 

EGRT 

rad 

DRAG 

Current drag acceleration 

State vector 

ft/sec 2 

EGFLG 

Guidance mode flag 

Entry processor 
executive 

n .d. 

HLS 

Altitude above runway 

State vector 

ft 

LOD 

Current lift/drag ratio 

Aerodynamics 

n.d . 

RDOT 

Current oblate Earth 
altitude rate 

State vector 

ft/sec 

ROLL 

Current bank angle 

Autopilot 

rad 

START 

Initialization flag 

Entry processor 
executive 

n.d. 

TRANGE 

Current range to runway 

EGRT 

n .mi . 

VE 

Current relative velocity 

State vector 

ft/sec 

VI 

Current inertial velocity 

State vector 

ft/sec 

XLFAC 

Current load factor 

Aerodynamics 

ft/sec 2 

MM304 

Preentry bank angle 

Med input 

deg 

MM3 Oil 

Preentry angle of attack 

Med input 

deg 

MEP 

Left-hand HAC select flag 

n .d . 
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TABLE 3-3-2.- ENTRY GUIDANCE OUTPUTS 


Symbol 

Description 

Unit 

Destination 

ALPCMD 

Angle-of-attack command 

rad 

Autopilot 

ROLLC 

Bank command 

rad 

Autopilot 

DREFP 

Drag reference 

ft/sec^ 

Display 

DRAG 

Actual drag 

ft/sec^ 

Display 

ROLREF 

Bank reference 

deg 

Display 

ISLECT 

Guidance phase indicator 
relative 

n.d. 

Display 

VCG 

Velocity at start of constant 
drag phase 

ft/sec 

Display 

VRR 

Velocity at first bank 
reversal 

ft/sec 

Display 

EOWD 

Energy overweight 

ft 

Display 

EEI 

Entry evaluation indicator 

n.d. 

Display 

RC176G 

First roll command after 0.176g 

deg 

Display 
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TABLE 3-3-3.- ENTRY GUIDANCE INTERNAL PARAMETER DEFINITIONS 


Symbol 

Description 


ACLAM 

Maximum allowable alpha 


ACLIM 

Minimum allowable alpha 


ACMD1 

Scheduled angle of attack 


ALDCO 

Temporary variable in phase 3 reference parameters 

ALDREF 

Vertical L/D reference 


ALPCMD 

Angle--of-attack command 


ALPDOT 

Rate of change of ALPCMD 


ARG(1 ) 

Cosine of commanded bank angle 


ARG(2) 

Cosine of unlimited bank command 


ARG(3) 

Cosine of bank reference angle 


A2 

Temporary variable used in computing 

range and updating D23 

CAG 

Pseudoenergy /mass used in transition 

(L/D) reference 

CQ1 (1,2) 

Constants in Ith temperature control 

D-V quadratic 

CQ2( 1,2) 

VE coefficients in Ith temperature control D-V quadratic 

CQ3( 1,2) 

p 

VE coefficients in Ith temperature 

control D-V quadratic 

Cl 

dD/dE in transition 


C16 

d(L/D)/dD 


Cl 7 

d( L/D)/dH 


C2 

Component of L/D reference 


cn 

Reference altitude rate term 


C20 

d ALPHA/dCd gain 


DD 

Drag - DREFP 


DDS 

Limited value of DD 
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TABLE 3.3-3.- ENTRY GUIDANCE INTERNAL PARAMETER DEFINITIONS - Continued 


Symbol 


Description 


DDP 

DELALF 

DELALP 

DLRDOT 

DLIM 

DLZRL 

DRDD 

DREF(1 ,2) 

DREFP 

DREFPT 

DREFP 1 

DREFP3 

DREFP4 

DREFP5 

DRF 

DX( 1,2) 

DZOLD 

DZSGN 

D231 

EEF 

HDTRF ( 1 ) 

I ALP 


Past value of DD 
Delta ALPHA from schedule 
Command ALPHA increment 
R feedback term 

Maximum value of DREFP in transition 
Test variable in bank-angle computation 
Derivative of range wrt drag 

DREFP for Ith temperature control D-V quadratic 

Drag reference used in controller 

DREFP-DF in transition phase 

DREFP in equilibrium glide 

DREFP test value for transition to phase 3 

DREFP test value for transition to phase 4 

DREFP test value for transition to phase 5 

Test value for transition to D23-VB1 quadratic reference 
parameters 

Normalized values of DREFP 
Previous value of DELAZ 
Change in DELAZ 
First updated value of D23 
Energy/mass 

Intermediate calculation of temperature control R ref 
ALPCMD segment counter 
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TABLE 3.3 
Symbol 
ICT 

ISLECP 

ISLECT 

ITRAN 

LMFLG 

LMN 

LODV 

LODX 

QO ,2,3) 

RCG 

RCG1 

RDEALF 

RDTREF 

RDTRF1 

REQ1 

RER1 

RF ( 1 ,2) 

RFF1 
RDTRF 
ROLLC(I) 
ROLLC(2 ) 
ROLLC(3) 

RPT 


-3.- ENTRY GUIDANCE INTERNAL PARAMETER DEFINITIONS - Continued 


Description 


Alpha modulation flag 

Past value of ISLECT 

Entry guidance subphase counter 

Transition initialization flag 

Saturated roll command flag 

Maximum value of LODV 

Vertical L/D command 

Unlimited vertical L/D command 

DREFP/VE in temperature control phase 

Constant drag phase range 

Constant component of RCG 

Roll bias for alpha modulation 

Altitude rate reference 

RDTREF in Phase 3 

Equilibrium glide range x D23 

Transition phase range 

Ith range segment in temperature control phase x D23 

Temperature control range x D23 

Altitude rate reference corrected for Cd 

Roll-angle command about body axis 

Unlimited roll command 

Rolref 

Desired range at VQ 


17 


80FM23 


TABLE 3.3-3.- ENTRY GUIDANCE INTERNAL PARAMETER DEFINITIONS - Concluded 


Symbol 


Description 


R231 

START 

T1 

T2 

T2D0T 
T20LD 
V(1 ,2,3) 

VB2 

VCG 

VE2 

VF(1 ,2) 

VSAT2 
VTRB 
VX(1 ,2) 

XLOD 

YL 

ZK 


Phase 2 and 3 range x D23 

Entry guidance first pass flag 

Equilibrium £vide vertical lift acceleration 

Constant drag level to reach target 

Rate of change of T2 

Old T2 value 

Velocity sampling points for temperature control numerical 
range prediction 

VB 2 

Phase 3-4 boundary velocity 
VE 2 

Upper velocity bounds for Ith temperature control 
range segment 

VSAT 2 

R feedback velocity lockout 

Velocities where dD/dV=0 in Ith temperature control D-V 
quadratic 

Limited iift/drag ratio 

Maximum heading error absolute value 

R feedback gain 
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TABLE 3.3-4.- AUTOPILOT INPUT DATA 
(a) Input parameters 


Symbol 

Description 

Source 

Unit 

ALPWWD 

Angle of attack with winds 

Integrator 

rad 

ALPHAP 

Previous pass angle of attack 

Entry integrator 

rad 

ALPCMD 

Angle-of-attack command 

Entry guidance 

rad 

PQR(IV1 ) 

Initial bank rate 

Entry integrator 

rad /sec 

PQR(IV+2) 

Initial alpha rate 

Entry integrator 

rad/ sec 

ROLLC 

Bank command 

Entry guidance 

rad 

ROLLP 

Previous pass bank angle 

Entry integrator 

rad 

RRPAST 

Previous pass bank rate 

Entry integrator 

rad/sec 

PRPAST 

Previous pass angle-of-attack 
rate 

Entry integrator 

rad /sec 

EGFLG 

Guidance mode flag 

Executive 

n.d. 

MM304 

Preentry bank command 

Med input 

rad 

MM304a 

Preentry angle-of-attack 
command 

Med input 

rad 

ROLL 

Current bank angle 

Integrator 

rad 

DELTT 

Integration time step 

Executive 

sec 
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TABLE 3.3-4.- AUTOPILOT INPUT DATA - Concluded 
(b) Input constants 


Symbol 

Description 

Value 

Unit 

Class 

RA 

Maximum bank acceleration 

0.029671 

p 

rad/ sec 

f 

RADB 

Bank attitude deadband 

0.00174533 

rad 

f 

RADB2 

Bank attitude deadband 
for EGFLG=2 

0.00174533 

rad 

f 

RRM 

Maximum bank rate 

0.087266 5 

rad/ sec 

f 

PA 

Maximum angle-of-attack 
acceleration 

0.0174533 

O 

rad/ sec 

f 

PADB 

Angle-of-attack attitude 
deadband 

0.00174533 

rad 

f 

PRM 

Maximum angle-of-attack 
rate 

0.0872665 

rad/sec 

f 

RAD160 

160 degrees converted to 
radians 

2.7925268 

rad 

f 

RAD180 

1 80 degrees converted to 
radians 

3.1415927 

rad 

f 

RAD360 

360 degrees converted to 
radians 

6.2831853 

rad 

f 

DTR 

Degrees to radians 

0.017453292 

rad 

f 

KONTROL 

Initialization flag 

15 

n. d. 

f 
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TABLE 3.3-5.- AUTOPILOT OUTPUTS 


Symbol 

Description 

Unit 

Destination 

ALPWND 

Updated angle of attack 
with winds 

rad 

Integrator 

BANK 

Undated bank angle 

rad 

In tegrator 

PR 

Updated angle-of-attack 
rate 

rad/ sec 

Integrator 

RR 

Updated bank rate 

rad/sec 

Integrator 
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TABLE 3.3-6.- AUTOPILOT INTERNAL PARAMETER DEFINITIONS 


Symbol 

Description 

BANK 

Bank about the velocity vector 

DIFF 

Difference between ROLLC and ROLL 

DIR 

Direction of acceleration 

DTIM 

DAP3D time step 

EFRATE 

Average value of EFRAT1 and EFRAT2 

EFRAT1 

DIFF rate 

EFRAT2 

Past value of EFRAT1 

ICPPLN 

Cycle frequency of DAP3D 

ICTTRN 

First and third pass flag 

INTRY 

Flag to establish constant set 

KONTRL 

PHSPLN initialization flag 

RR1 

Current bank rate 

ROLL 

Current bank angle 
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TABLE 3.3-7.- TARGETING ROUTINE INPUT DATA 
(a) Input parameters 


Symbol 

Description 

Source 


Unit 

AZRW 

Bearing from true north of 
runway +X axis 

Landing 

site 

table 

rad 

ALTD 

Altitude of runway 

Landing 

site 

table 

ft 

RLAT 

Redesignation runway latitude 

Landing 

site 

table 

rad 

RLONG 

Redesignation runway longitude 

Landing 

site 

table 

rad 

RALTD 

Redesignation runway altitude 

Landing 

site 

table 

ft 

RAZRW 

Redesignation runway azimuth 

Landing 

site 

table 

rad 

TLONG 

Runway longitude 

Landing 

site 

table 

rad 

TLATD 

Runway latitude 

Landing 

site 

table 

rad 

VREDS 

Redesignation velocity 

Landing 

site 

table 

ft/sec 

XYZE 

Vehicle position vector 

Integrator 


ft 

XYZED 

Vehicle velocity vector 

Integrator 


ft/sec 

SLATD 

Secondary runway latitude 

Landing 

site 

table 

rad 

SLONG 

Secondary runway longitude 

Landing 

site 

table 

rad 

SRTE1 

Secondary runway altitude 

Landing 

site 

table 

ft 

SRAZ 

Secondary runway azimuth 

Landing 

site 

table 

rad 
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TABLE 3.3-7.- TARGETING ROUTINE INPUT DATA - Concluded 


(b) Input constants 


Symbol 

Description 

Value 

Units 

Class 

RTURN 

Radius of heading alinement 
circle 

20 000 

ft 

m 

XNEP 

Distance between HAC center 
and threshold point 

-39 519 

ft 

m 

RX22 

p 

Polar radius /equatorial 
radius 2 

0.9933065782 

n.d 

f 

VMIDPT 

Velocity limit for midpoint 
targeting 

30 000. 

fps 
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TABLE 3-3-8.- TARGETING ROUTINE OUTPUTS 


Symbol 

Description 

Unit 

Destination 

TRANGE 

Range to target 

n. mi. 

Entry guidance 

DELAZ 

Azimuth error 

rad 

Entry guidance 

RCHMAG 

Radius of landing site 

ft 

Deorbit processor 
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TABLE 3-3-9.- TARGETING ROUTINE INTERNAL PARAMETER DEFINITIONS 


Symbol 

Description 


RLS 

Landing site in Earth-fixed coordinates 


RC 

Center of heading alinement circle in runway coordinates 

HACEF 

Center of heading alinement circle in Earth-fixed 

coordinates 

BARCC 

Heading to center of alinement circle 


CTHVC 

Cosine (angle between vehicle and HAC) 


BARWPI 

Heading to tangent point on HAC 


DARC 

Distance around heading alinement circle 


TT2DIN 

Earth- fixed to top ode tic matrix 


PSI 

Heading of vehicle 


REC 

Runway to Earth- fixed matrix 


IP ASS 1 

Midpoint targeting flag to denote computation of 
runway range and delta azimuth 

primary 

XLAT 

Primary runway latitude save parameter 


XLONG 

Primary runway longitude save parameter 


XRTE1 

Primary runway altitude save parameter 


XRAZ 

Primary runway azimuth save parameter 
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3.4 ENTRY GUIDANCE FORMULATION 


3.4.1 Requirements Overview 

The entry guidance is the source for bank-angle and angle-of-attack commands, 
which are used to control the entry trajectory. The entry guidance can be 
called by either of two modes. If the entry guidance flag (EGFLG) is set to 
zero or one, a normal entry guidance function is exercised that will duplicate 
the Orbiter avionics entry guidance function. This mode will be used to simu- 
late an entry trajectory and will also be used in the final iteration mode of 
the entry target generator (ETG) processor. The second mode (EGFLGr2) is used 
in the entry processor to simulate an entry based on a "canned" drag profile for 
the ETG targeting processor. 

In the normal mode (EGFLG=0), the entry guidance controls the entry trajectory 
by bank-angle modulation while using a preselected angle-of-attack profile, 
which is a function of relative velocity. Range predictions are based on solu- 
tions to the equation of motion for a specified entry drag-velocity profile. The 
drag-velocity profile, the shape of which is specified by the mission cons ants 
table, consists of quadratic, pseudoequilibriura glide, linear, and constant drag 
segments. Downrange errors are nulled by changing the magnitude of the bank 
angle, and crossrange errors are nulled by bank-angle reversals, which limit the 
crossrange error within a converging deadband. This mode begins at 400 000 ft 
and ends at TAEM interface, relative velocity = VTAEM. 

In the ETB "canned" mode (EGFLG=2), the entry guidance controls the entry 
trajectory to a predefined drag-velocity profile. The entry simulations begin 
at 400 000 ft and are terminated after the pullout maneuver has been completed 
(and the flight profile stabilized) on the drag velocity profile at exactly 23 
000 ft/ sec relative velocity (VETG) . 

In the normal mode, the entry guidance consists of five major phases: the 

preentry phase, the temperature control phase, the equilibrium glide phase, the 
constant drag phase, and the transition phase, as shown in figure 3 -4. 1-1. 

The preentry phase maintains the vehicle attitude until a predetermined total 
load factor level is reached (ASTART). During this phase, the vehicle is 
maintained in a three-axis attitude hold mode. This attitude is specified by 
the flight controller as MM304<J>, MM304a, and MM304g. This phase is terminated 
at 5.66 ft/sec^ total acceleration level and the temperature control phase is 
entered . 

The temperature control phase is designed to control the entry trajectory 
through pullout to a temperature profile consistent with the overall entry pro- 
file shape optimization and entry ranging requirements. This phase consists of 
two quadratic drag-velocity segments, which are used to optimize the entry pro- 
file. Range predictions are based on analytic solutions of the equations of mo- 
tion for these two segments. This phase is terminated at the velocity, VB1 
(specified in the mission constants table), and the guidance is transferred to 
the equilibrium glide phase at this point. 
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Hie equilibrium glide phase produces an equilibrium glide-type trajectory, con- 
sistent with the ranging solution, until the equilibrium glide trajectory inter- 
sects the constant drag trajectory required to reach the target. At this point 
control is transferred to the constant drag phase, and range predictions are 
based on a constant drag profile until the transition phase is entered. At a 
specified velocity, VTRAN (mission constants table), control is transferred to 
the transition phase. The transition phase is based on a linear drag profile, 
as a function of energy, which is required to null the range error. The transi- 
tion from the entry guidance to the TAEM guidance occurs at a velocity, VTAEM, 
specified in the mission constants table. Presently, this transfer point is at 
an Earth-relative velocity of 2500 ft/sec. 

'Ihe entry guidance generates bank-angle and angle-of-attack commands to be used 
by the autopilot . Ihe bank-angle commands are designed to converge the actual 
drag acceleration level to the reference drag-velocity profile, described above, 
that is consistent with the ranging solution. The bank-angle command is 
generated from a vertical L/D command (which is a function of a reference L/D) , 
the difference between drag and drag reference, and the difference between alti- 
tude rate and altitude-rate reference. The angle-of-attack profile is a func- 
tion of Earth-relative speed and consists of a series of linear and quadratic 
segments. The angle-of-attack profile is controlled through inputs in the mis- 
sion constants table. A complete derivation of all entry guidance equations can 
be found in reference 3- 


< 
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3-4.2 Entry Guidance Executive (EGEXEC) 

The EGEXEC function calls the other entry guidance functions in proper sequence 
and controls guidance phase transitions. In the normal mode (EGFLG=0 or 1 ) , 
entry guidance is divided into five phases: preentry, temperature control, equi- 

librium glide, constant drag, and transition. Each phase, except preentry, com- 
putes a reference drag acceleration profile, which is based upon the ranging re- 
quirement and the vehicle constraints. 


The computations performed during each guidance cycle vary with the guidance 
phase . The active guidance phase is defined by the integer variable ISLECT ac- 
cording to the following table. 


ISLECT 

1 

2 

3 

4 

5 


Guidance Phase 
Preentry 

Temperature control 
Equilibrium glide 
Constant drag 
Transition 


The entry guidance phase transfer logic is summarized in table 3.4.2. 
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'Ihe preentry phase commands a preselected bank angle in the EGPEP function. The 
preentry phase (ISELECT = 1) can be terminated by any of three conditions. 

a. Normal termination: The preentry phase is normally terminated at 5.66 

ft/sec 2 by setting the ISLECT flag to 2 when the total acceleration exceeds 
the threshold value of ASTART ft/sec 2 . 

b. Alternate termination: If, at the threshold load factor level, the current 

constant drag level to reach the target is greater than the desired constant 
drag level (ALFM), ISLECT is set to 4. This could occur for an extremely 
short-range case. 

c. Threshold load factor level: If the current relative velocity is less than 

VTRAN , ISLECT is set to 5. 

The temperature control phase (ISLECT = 2) computes the required drag velocity 
profile during the high-temperature region of entry. The functions EGRP and 
EGREF contain the detailed temperature control equations used by the entry guid- 
ance . During the temperature control phase , the drag-velocity reference 
trajectory is divided into two quadratic segments; function EGREF determines 
which of these quadratics is to be used. Figure 3. 4. 2-2 presents example qua- 
dratics for this phase. The quadratic for the higher velocity region is used 
when 


VE > VA1 


and 


DRF = (DREF(2) - DREF(l)) (DREF(2) - DREF (1)) + (HDTRF(I) 
- HDTRF(2) GS1 ) > 0 


where DREF ( 1 ) , HDTRF(I), DREF (2) , and HDTRF(2) are all computed in EGREF. If 
either test is failed, then the reference parameters RDTREF and DREFP 
computed for the lower velocity drag profile are used. The quadratic seg- 
ment switching point is controlled by the mission load parameter, GS1 . The 
quadratic switching logic, in the EGREF function, is similar to that used 
for guidance phase transitions. 

The temperature control phase can be terminated in four ways. 

a. Normal termination: The temperature control phase normally transfers to 

the equilibrium glide phase, ISLECT = 3. Transfer is planned to occur when 
the drag reference profiles for the temperature control phase and the equi- 
librium glide phase intersect, as illustrated in figure 3 .4. 2-3. However, 
if the slopes of the two intersecting profiles are different, the bank-angle 
command may jump at the transfer point . Bank-smoothing logic provides a 
smooth roll-angle command history for minimum reaction control system (RCS) 
fuel usage by transferring phases when the commanded bank angle is the same 


31 



Figure 3. 4. 2-2.- Temperature control phase quadratic definition. 



Figure 3. 4. 2-3.- Bank-angle smoothing logic between guidance phases 
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■ for both phases. This transition occurs before the intersection point 
of the drag profiles and is accomplished using the following equation: 


If VE < VA and 

DREFP < DREFP3 

then 

ISLECT = 3 


where DREFP3 = DREFP1 + GS2 (RDTREF - RDTRF1) and DREFP3 is computed in the 
EGREF function. The value of GS2 determines the transfer point. 

b. Alternate termination: The guidance transfers to phase 3 when VE < VB1 if 

the bank-smoothing logic has not been satisfied by that time. Thus, if VE 
< VB1 , then ISLECT = 3. 

c. Alternate termination for short ranges: For very short-range targets, the 
desired constant drag level may be reached before equilibrium glide phase 
initiation. In this case, control is transferred directly from the tempera- 
ture control phase to the constant-drag phase (ISLECT = 4) when 


VE < VCG + DELV 


and 


DREFP > DREFP 4 

where DREFP4 = T2 + GS3 (RDTREF + 2 HS T2/VE) and DREFP4 is computed in 
EGREF. Bank-command smoothing is provided through the constant GS3. 

d. Extremely short-range termination: Transition to phase 4 occurs 

If T2 > ALFM 

then 


ISLECT = 4 

where T2 is computed in EGCOMN. 

The equilibrium glide phase (ISLECT = 3) shapes the drag velocity profile so 
that the constant drag level to reach the target converges on the desired 
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constant drag level, ALFM. The functions EGRP and EGREF contain the equilibrium 
glide equations. There are three possible transfers from the equilibi ium glide 
phase . 

a. Normal termination: The equilibrium glide phase transfers to the constant 

drag phase, ISLECT = 4. Bank-smoothing logic is provided by the constant 
GS3. Transfer occurs when 


VE < VCG + DELV 


and 


DREFP > DREFP4 


where DREFP4 = T2 + GS3 (RDTREF + 2 HS T2/VE) and is again computed in 
EGREF. 

b. Alternate termination: For very long-range trajectories, the predicted ve- 

locity at the intersection of the equilibrium glide and constai t drag phases 
is less than the transition phase initiation velocity, VTRAN. When this 
occurs, the equilibrium glide phase transfers directly into the transition 
phase, ISLECT = 5. The transition to phase 5 occurs when 


VE < VTRAN + DELV 


and 


VCG < VTRAN 


and 


DREFP > DREFP5 


where DREFP5 = DF + (EFF - EEF4) Cl + C-S4 (RDTREF - RDTRFT) is computed 
in EGCOMN. The variables EEF, Cl, and RDTRFT are also computed in EGCOMN. 
As in other phases, bank-smoothing logic is provided by the mission load 
constant, GS4. 

c. Alternate termination: For extremely short ranges, as in the temperature 

control phase, transition to the constant drag phase (ISLECT = 4) occurs 
if 
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T2 > ALFM 


where T2 is again computed in EGCOMN. 

The constant drag phase (ISLECT = 4) shapes the entry profile along a constant 
drag velocity profile to maximize the control system margins. Function EGREF4 
contains the constant drag range prediction and controller equations. The con- 
stant drag phase terminates and transfers to the transition phase (ISLECT = 5) 
at a predetermined time before the transition energy level, ETRAN. Transfer 
occurs when 


VE < VTRAN + DELV 


and 


DREFP > DREFPS 


where DREFP5 = DREFP5 in the equilibrium glide phase and is computed in EGCOMN . 

The transition phase is the final entry phase and is used to 3teer the Orbiter 
to the proper TAEM interface conditions. The transition phase and entry guid- 
ance are terminated at the start of the TAEM major mode. The entry-to-TAEM tran- 
sition logic is defined by 


If (ISLECT /I and VE < V TAEM) then EG END = 1 


To execute the entry guidance computations properly, the following functions 
must be called in the sequence shown. 

a. The EGSCALEHT function is called. 

b. On the first entry guidance pass (START = 0), EGINIT is then executed. 

c. The EGCOMN function is called. 

d. The phase transition logic for the preentry phase (ISLECT = 1) and the first 
alternate termination test for the temperature control phase (if VE < VB1 , 
then ISLECT = 3) are executed within EGEXEC. The tests to transfer to phase 
4 if T2 > ALFM are also made at this time. 
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The functions next called depend on the value of ISLECT. 


If ISLECT = 
If ISLECT = 
If ISLECT = 
If ISLECT = 


1 , then perform EGPEP 

2 or 3, then perform EGRP, and then EGREF 

4, then perform EGREF4 

5, then perform EGTRAN 


After these functions have been executed, the remainder of the phase transition 
logic may be performed with EGEXEC at any time. The output commands are then 
computed by calling EGALPCMD (angle-of-attack command) and the sequence 
EGGNSLCT , EGLODVCMD , and EGROLCMD (bank-angle command). If ISLECT = 1, EGGNSLCT 
and EGLODVCMD are bypassed, and only EGALPCMD and EGROLCMD are called. For all 
values of ISLECT, either the angle-of-attack or bank-angle command may be 
computed first. The entry guidance execution sequence is summarized in figure 
3. 4. 2-4. 

In the ETG "canned" mode (EGFLG = 2), ISLECT is never allowed to be greater than 

2 . 


3.4.3 Entry Guidance Scale Height (EGSCALEHT) 

The guidance function EGSCALEHT generates an altitude scale height (of 
atmospheric density) modeled on the 1962 standard atmosphere. This parameter is 
used in calculating the altitude rate reference term 


RDTRF = -HS — + - + — 

VE D cd 

where D = drag acceleration 

D = time derivative of D 
cd = drag coefficient 
cd = time derivative of cd 

Empirical curve fits of the altitude scale height (HS) as a function of relative 
velocity (Ve) have been implemented into the entry guidance. 
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Figure 3.^.2-^.- Entry guidance 

sequence . 
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3.4.4 Entry Guidance Initialization (EGINIT) 

The guidance function EGINIT serves as the initialization routine for entry guid- 
ance. In this routine initial values are set, and parameters calculated only 
one time are computed. 


3.4.5 Entry Guidance Common Computation (EGCOMN) 

The entry guidance contains several parameters used continuously throughout the 
guidance program. These are computed in EGCOMN, and are such parameters as 
energy (EEF) the constant drag level to reach the target (T2), and the rate of 
change of 12 ( '12 DOT). 


3*4.6 Entry Guidance Preentry Phase (EGPEP) 

In the Orbiter avionics system, the purpose of EGPEP is to generate a vertical 
L/D command (LODV) by means of the ILOAD parameter PREBNK. However, in order 
for the MCC to simulate either an automatic or manual preentry phase, the LODV 
equation in the Orbiter avionics system should be replaced by the bank-angle 
input in MM3Q4<j> by the flight controller. Also, the angle-of-attack command 
issued in EGALPCMD should be overridden by MM304a. The preentry phase (ISLECT 
= 1) is terminated by EGEXEC at a sensed total acceleration level equal to 
AS TART (currently 5.66 ft/sec 2 ). 


3.4.7 Ehtry Guidance Range Prediction (EGRP) 

The EGRP function serves as the range predictor during the temperature control 
and equilibrium glide phases. The range prediction is then used to determine 
the proper drag level during phases 2 and 3 to achieve the desired range at the 
entry- TAEM interface and is only called when ISLECT is equal to 2 or 3. 

In order to determine the proper drag-velocity profile for range control, a 
range prediction is made of the entire entry trajectory. This is accomplished 
by using five drag-velocity segments; two during the temperature control phase, 
one during the equilibrium glide phase, one during the constant drag phase, and 
a constant range value during the transition phase. The temperature control, 
equilibrium glide, and constant drag phase range segments are computed in EGRP. 
The range value for transition, RPT, is computed in EGINIT and is constant in 
order to provide a nominal transition range at the transition interface, VTRAN. 
Once the transition phase is entered, the transition range prediction is 
modified in the EGTRAN function to meet the range requirements. 

The drag-velocity segments during the temperature control, equilibrium glide, 
and constant drag phases are anchored at specific velocity points, as illustrated 
in figure 3 -4 .7-1 • 

The range for the temperature control phase is predicted along two quadratic 
drag-velocity segments anchored at VB1 and at VA1. For the equilibrium glide 
segment, the range is predicted between VB1 and VCG, the computed intersection 
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Figure 3. 4. 7-1.- Drag-velocity segments used in range predictions 
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Figure 3. 4. 7-2.- Freezing of equilibrium glide profile, 
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point between the equilibrium glide phase and constant drag phase, and for the 
constant drag segment, the range is predicted between VCG and VQ. In all cases, 
the entire drag-velocity profile is anchored at a drag level of D23 at a veloc- 
ity of VB. During the temperature control phase, VB is defined as VB1 , and dur- 
ing the equilibrium glide phase, VB is defined as the current relative velocity. 
Therefore , 


If 

VE 

> VB1 

VB = 

VB1 

If 

VE 

< VB1 

VB = 

VE 


In the equilibrium glide phase, as the drag profile approaches the desired con- 
stant drag level, the focus of the equilibrium glide drag reference parameter 
may wander from a precise equilibrium glide profile shape trying to drive T2 to 
precisely ALFM, the desired constant drag level. In order to provide a more 
uniform equilibrium glide drag profile at the junction point between the equilib- 
rium glide and the constant drag phases, the equilibrium glide reference profile 
is "frozen" when the rate of change of T2 is near zero, as illustrated in figure 
3. 4. 7-2. This is accomplished by freezing VB at the current value of VE when 
T2D0T becomes less than DT2MIN and when VE is less than VCG + DELV. 

If EGFLG is equal to 2 (ETG canned mode), the ranging iteration of D23 is 
bypassed and D23 is set equal to D23C in order to provide a "canned" drag veloc- 
ity profile for the ETG mode. 


3.4.8 Entry Guidance Reference Parameters (EGREF) 

In order to control the Orbiter to the desired drag-velocity profile required 
for range control, a bank-angle command is generated from a vertical L/D command 
equation. This vertical L/D command equation consists of actual and reference 
parameters. The function of EGREF is to generate the drag reference, the alti- 
tude rate reference, and the phase-dependent part of the L/D reference parameter 
for the temperature control and the equilibrium glide phase. This function is 
only called when ISLECT is equal to 2 or 3. 


3.4.9 Entry Guidance Constant Drag Phase (EGREF4) 

The purpose of the EGREF4 function is to generate the drag reference, the alti- 
tude rate reference, and the phase-dependent part of the L/D reference for the 
constant drag phase. This function is called only when ISLECT is equal to 4. 


3.4.10 Entry Guidance Transition Phase (EGTRAN) 

The transition phase function (EGTRAN) computes the range potential from the drag 
reference level at the end of the constant drag phase to the transition phase 
end target conditions (DF and EEF4) and then computes the correct drag energy 
profile to null any range error. EGTRAN also computes the controller reference 
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parameters : drag reference , altitude rate reference , and the phase-dependent 

part of L/D reference. This function is called only when ISLECT is equal to 5. 

3.4.11 Entry Guidance Angle-of-Attack Function (EGALPCMD) 

The EGALPCMD function generates the angle-of-attack command for the flight con- 
trol system. The angle-of-attack profile commanded by the entry guidance is a 
preselected profile established by means of mission-dependent constants. The 
entry velocity regime is divided into NALP+1 segments, and the commanded angle 
of attack in each segment is defined by a quadratic function of relative veloc- 
ity. Figure 3.4.11-1 shows a typical angle-of-attack command profile and illus- 
trates the flexibility available in the profile selection. 


3.4.12 Entry Guidance Gain Selection Function (EGGNSLCT) 

The EGGNSLCT function computes the drag error gain (C16) and the altitude rate 
error gain (C17) in the controller vertical L/D command equation. 'Ihese gains 
are a function of the actual drag acceleration level and the difference between 
drag and drag reference. 


3 . 4.13 Entry Guidance Lateral Logic and Vertical L/D Command Function 

( EGLODVCMD) 

The purpose of the EGLODVCMD function is to 

a. Compute the L/D reference parameter (ALDREF) . 

b. Compute an R feedback term to correct drag error biases caused by poor navi- 
gation (DLRDOT). 

c . Compute the vertical L/D command from the controller equation 

LODX = ALDREF + Cl6 (DRAG - DREFP) + C17 ( RDTRF + DLRDOT - RDOT) 

d. Perform a velocity check to see if angle-of-attack modulation should begin 
in order to keep drag on the drag reference profile. 

e. Compute the bank angle limit (LMN) and, finally, to compute the bank direc- 
tion (RK2R0L) . 


3.4.14 Entry Guidance Bank Command Function (EGROLCMD) 

The purpose of the EGROLCMD function is to generate a bank command for the 
autopilot and a bank-reference parameter for display. The bank command is 
computed from the vertical L/D command parameter and the bank reference is 
computed from L/D reference. If the angle-of-attack modulation flag (ICT) is 
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equal to one, a bank-command bias is computed as a function of the ALPHA 
difference with respect to the nominal ALPHA schedule. 


3-4-15 Eh try Guidance Data Flow Summary 

The data flow charts (figs. 3.4.1 5—1 and 3.4.1 5—2) present the data flow of all 
computed and stored parameters within the entry guidance. 


3.5 ENTRY AUTOPILOT FORMULATION 


3.5.1 Requirements Overview 

A simple three-degree-of-freedom autopilot is required to execute the bank- 
angle and angle-of-attack commands generated by the entry guidance. The 
autopilot requirements are the same for both of the entry guidance modes. 

The autopilot, by means of a simple phase plane, generates the actual bank angle 
and angle of attack based on a maximum acceleration and rate limit in the bank 
and pitch channels. The sideslip angle (8) is always assumed to be zero. 

The entry guidance generates a bank-angle and angle-of-attack command necessary 
to control the trajectory. The autopilot generates the Orbiter attitude re- 
sponse to these commands over the next computer cycle (2.0 sec), ignoring the 
high-frequency dynamics. Eased on the attitude response characteristics, the 
autopilot determines if the commanded attitude can be achieved during the next 
computer cycle, and if commanded attitudes cannot be achieved the autopilot de- 
termines the achievable attitude at the end of the computer cycle. If the atti- 
tude can be reached within the computer cycle, a deadband attitude and rate is 
established about the commanded attitude. This new attitude is then used to com 
pute the trajectory dynamics and the accelerations for the next integration step 
during entry. 


3.5.2 Autopilot Executive (DAP3D) 

The autopilot executive routine (DAP3D) is the driver routine for the simplified 
autopilot phase plane (PHSPLN). Assuming a four-pass Runge-Kutta integrator, 
PHSPLN is called on the first and third pass by means of the pass counter 
(ICTTRN). ROLPLN and PCHPLN are entry points for the generalized routine PHSPLN 
and are called for bank-attitude control and angle-of-attack control, respec- 
tively. 

DAP3D also computes the Orbiter attitude with respect to the velocity vector, 
bank, sideslip, and angle of attack. These attitudes are then used in the gener 
ation of the body-to-inertial coordinate system transformation matrix. Table 
3.3-4 presents the inputs to the autopilot, table 3.3-5 presents the outputs, 
and table 3.3-6 presents the internal parameter definitions. Appendix B pres- 
ents the formulation flow charts, and appendix D presents the IBM structured 
flow charts for the autopilot. 
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Figure 3. 4.15-1.- Entry guidance external data flow summary. 
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Figure 3.4.15-2.- Entry guidance internal data flow. 
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Figure 3.4.15-2.- Continued. 
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STORAGE 
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ALPCMD, ALPDOT 
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| DLRDOT, ALPHA, 

! ALDREF, LODX, EEI 
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Figure 3. A. 15-2.- Continued. 
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Figure 3. 4. 1 5—2. — Concluded. 
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3.5.3 Autopilot Phase Plane (PHSPLN) 

The PHSPLN routine is a simplified phase plane that is used to generate the 
Orbiter bank-angle and angle-of-attack attitude. The routine PHSPLN has two 
entry points (ROLPLN and PCHPLN), which update the bank-angle and angle-of- 
attack attitudes, respectively. The phase plane logic uses an acceleration 
level (RA or PA), a maximum rate (RRM or PRM) , an attitude deadband (RADB or 
PADB), and a bank-angle or angle-of-attack error to determine the current 
Orbiter attitude. A variable (ICPPLN) has been added to allow 'he calling fre- 
quency of the routine, per integration At, to be selected by . ' i user. 


3.5.4 Autopilot Data Flow Summary 

The data flow charts (figs. 3 .5. 4-1 and 3. 5. 4-2) present the data flow of all 
computed and stored parameters with the autopilot. 


3.6 TARGETING ROUTINE (EGRT) 

The targeting routine, EGRT, computes the great circle range from the Orbiter to 
the runway threshold point via the heading alinement circle. This is accom- 
plished by determining the tangent point or. the heading alinement circle of r 
vector from the vehicle to the alinement circle. This tangent point is 
converted into ar. Earth-fixed position vector, and the great circle range to tar- 
get is computed between the vehicle and this tangent point. The arc length is 
then computed from the tangent point around the alinement circle to WP1, the 
straight-in approach point. The range to target is then computed as the sum of 
the great circle range to the tangent point, the arc length around the alinement 
circle, and the distance between WP1 and the runway threshold point. The azi- 
muth error is computed as the difference between the vehicle Earth-relative azi- 
muth and the heading to the heading alinement tangent point. The targeting will 
nominally select the nearest heading alinement circle on the same side of the 
runway centerline as the Shuttle position. The option to force the selection of 
the left-hand heading alinement circle independent of Shuttle position is 
exercised by setting the HAC select flag (MEP=1) by MED input. Midpoint 
targeting is needed to minimize the delta range and azimuth error for a redesignation 
due to low altitude winds. Midpoint targeting i3 accomplished by targeting 
to both primary and secondary runways, and averaging the range and azimuth 
error values, and is terminated when 

a. The crew selects the desired runway 

b. The relative velocity satisfies a preset velocity limit and the primary 
runway is selected by default 

The flow charts for the targeting routine are found in appendix C and are 
subdivided into the following functions: 
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Figure 3-5. 4-1.- Autopilot external data flow. 
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Figure 3. 5. 4-2. - Autopilot internal data flow. 
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3.6.1 EGRT-EXEC, Targeting Executive Logic 

This routine calculates the Earth-fixed to runway transformation matrix for the 
primary and secondary runway (if performing midpoint targeting), sets the entry 
guidance redesignation and midpoint targeting flags, averages range and azimuth 
error for midpoint targeting, and calls the following subfunctions: 


3.6.2 EGRT-CHACRC, Center of Heading, Alinement Circle - Runway Coordinates 

This routine computes the position of the center of the heading alinement circle 
in runway coordinates. 


3 . 6.3 EGRT-CHACRC, Center of Heading, Alinement Circle in Earth-fixed 
Coordinates 


This function transforms the center of the heading alinement circle vector to 
the Earth-fixed coordinate system. 


3.6.4 EGRT-BV, Bearing of the Vehicle 

This function computes the heading of the vehicle based on the current vehicle 
Earth-fixed position vector. 


3.6.5 EGRT-BVCHAC, Bearing to Center of the Alinement 
Circle 


This function computes the heading of the vehicle to the center of the heading 
alinement circle. 


3 . 6.6 EGRT-COSTHETA, Great Circle Arc 

This function computes the great circle arc between the vehicle and the center 
of the heading alinement circle. 


3.6.7 EGRT-DWP 1 , Distance to WPI 

This function computes the range to the tangency point on the heading alinement 
circle. 


3.6.8 EGRT-DVNEP, Range- to-Threshold Point 

This function computes the heading to the tangency point on the heading 
alinement circle, the distance around the alinement circle, and the final range 
to runway threshold point. This function also saves the range to the primary 
runway in the case of midpoint targeting. 
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3.6.9 EGRT-DELAZ, Azimuth Error 

This function computes the azimuth error between the vehicle heading and the 
heading to the tangency point of the alinement circle. It also saves the azi- 
muth error to the primary runway in the case of midpoint targeting. 
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APPENDIX A 

ENTRY GUIDANCE FLOW CHARTS 

The following flow charts define the entry guidance formulations. 
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Fig ire A-1 . - Continued 
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Figure A-2. - EGSCALHT, scale height 
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INITIALIZE VARIABLES AND FLAGS 


CZOLD r 0.; 

I VRR s 0 
ITRAN = OFF; 

ISLECT = 1; 

ICT = 0 
IDBCHG = 0 
T2 = 0.; 

DREFP =0.; 

VQ2 = VQ VQ; 

RK2R0L = -SIGN(DELAZ); 

DLRDOT = 0. ; 

LMFLG = 0; 

VTRB = 60000.; 

DDP = 0.; 

RK2RLP = RK2ROL; 


! COMPUTE DESIRED TRANSITION RANGE ( RPT ) 

I 

RPT = -((ETRAN - EEF4) LOG(DF / ALFM) / ( 

ALFM - DF) + (VTRAN VTRAN - VQ2) / ( 

2. ALFM)) CNMFS + RPT 1; 

1 r " 

i 


VSAT2 = VSAT VSAT; 
VSIT2 = VS1 VS 1 ; 
VCG = VQ; 

D23 = D230; 

DX -| = 1 . ; 


,1: 

\ / 

Figure A-3.- EGINIT, initialization. Page 1 of 2 
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V0 1 = VB1; ! 

! 

vx 1 = VA; ! 

i 

VF 1 = VA1; i 

i 

A-, = AK; 1 

i 

V0 2 = VA 1 ; i 

i 

VX 2 = VA2; ! 

i 

A 2 = AK1; i 

i 

i 

VB 2 = VB 1 VB 1 ; I 

i 

i 

| COMPUTE C Of-! PON ENT OF CONSTANT DRAG PHASE RANGE, PCG 


i 

! RCG 1 = CNMFS (VSIT2 - VQ2) / (2. ALFM) ; 


IALP = NALP; 
START = 1; 


Figure A-3. - Concluded 
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COMPUTE COfflON VARIABLES 
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Figure A-4.- EGCOMN, comnon 
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COMPUTE VERTICAL L/D DURING PREENTRY PHASE 


LODX XLOD COS (MM30H* / RADEG) 
LODV = LODX; 


Figure A-5.- EGPEP, preentry phase. 
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Figure A— 6,— Concluded* 
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EGREF, reference parameters. 



Figure A-7-- Concluded. 
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COMPUTE REFERENCE PARAMETERS DURING CONSTANT DRAG PHASE 


! DREFP s T2; 


| RDTREF = -2. HS T2 / VE ; 


! DRDD = (TRANGE - RPT) / T2; 


i C2 = 0. ; 


i IT RAN = ON; 



Figure A-8.- EGREF4, constant drag phase. 
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COMPUTE REFERENCE PARAMETERS DURING TRANSITION PHASE 
THEN 


IF ITRAN = OFF 


\ 


/ 


DREFPT = DREFP - IF; 


DREFP = ALFM; i 
ITRAN = ON; 


\ THEN ! 

IF ABS (DREFPT) < El > i DREFP = DF + El SIGN (DREFPT ); 

/ I 


\ THEN ! 

IF DREFP < El > 1 DREFP = El; 

/ i 


DREFPT = DREFP - DF ; 

Cl = DREFPT / (EEF - EEF+); 

RER 1 = CNMFS LOG (DREFP / DF) / Cl; 

DRDD = MIN (SCALAR (CNMFS / (Cl DREFP) - 
RER 1 / DREFPT, DR DDL ) ) ; 

DREFP = DREFP + (TRANGE - RER 1 - RPT 1 ) / 
DRDD; 

DLIM = ALIM DRAG / XLFAC ; 

“ "™ — — 
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P 
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Figure A-9.-EGTRAN, transition phase. 
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\ THEN 

IF DREFP > DLIM > 

/ 


\ THEN 

IF DREFP < El > 

/ 


RDTREF = -AS VE (2. DREFP - Cl VE2) / 
CAG; 

C2 = 4. VE DREFP (1. - (VE2 /CAG) 2 ) / 
CAG; 


DREFP s El; 


DREFP = DLIM; | 

| 

Cl r 0.0; i 


Figure A-9.- Concluded. 
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COMPUTE ANGLE OF ATTACK COMMAND 


\ 

DO WHILE VE < VALP IALp AND IALP > 0; > 

/ 


IALP = IALP - 1; 


J = IALP + 1; 

ALPCMD = CALPOj + VE (CALPIj + CALP2jVE); 



“\ 

THEN 

1 I 

1 | 

! IF ISLECT = 1 
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1 ALPCMD = MM 304a | 
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t 1 

1 1 


ALPDOT = (ALPCMD - ACMD1 ) /DTEGD 
ACMD1 = ALPCMD 


Figure A-10.- EGALPCMD, angle-of-attack command. Page 1 of 1 
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COMPUTE CONTROLLER GAINS 


Cl 6 = CT16-) DRAG CT16 2 


\ THEN 

IF VE < VC16 >- 
/ 



C16 = MIDVAL (C16, CT16MN, CT16MX) 


IF ICT 



Cl 7 = MIDVAL (CT17 1 DRAG CT1 72 t CT17MN, 
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IF ICT = 1 > 

/ 



Figure A-11.- EGGNSLCT, gain select, 
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EGLODVCMD, lateral logic and vertical L/D command 
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Figure A-12.- Concluded. 
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APPENDIX B 

ENTRY AUTOPILOT FLOW CHARTS 

The Following flow charts define the entry autopilot formulations. 
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Figure B-2.- Concluded. 
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APPENDIX C 

TARGETING FLOW CHARTS 

The following flow charts define the targeting function for the entry processor. 
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Figure 0-1.- EGfiT-E/EC, targeting executive logic. Page 
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EGRT CHACRC 


RC(1 ) = XNEP ! 
RC( 3 ) = 0. i 


\ 



DO FOR I = 1,3 > 

/ 





[RG] = [REC][RGEF] 

RC( 2 ) = SIGN (RTURN, RG(2), 





\ 


IF MEP = 1 

> 



/ 


I 

i 

! 

i 


RGEF(I) = XYZE(I) - RLS(I) 


RC(2 ) = -RETURN 


■i ->t 


j: 


T) * 


Hh .h qij Ality 


Figure C-2.- EGRT-CHACRC, center heading alinement circle - 
runway coordinate system. 
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EGRT-CHACEFC 


[HACEF ] = [REC]“1[RC] 


i 




Figure C-3.- EGRT-CHACEFC, center heading alinement circle - 
Earth-fixed cooredinates. 
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EGRT-BV 


! \ 

1 DO FOR T - 1 R \ - , 

t 1 

I 1 

} RCOFFrn - RTS CM -t- HAfFFf'T'l ! 

1 LJ v* * UA X •• ’ | J ' 

i / 

t i 

i i 


I 

I 

I 

I 

I 

I 

\ 


VNORM(I) = XYZE(2)RCCEF(3) - XYZE (3 )RCCEF(2 ) 
VN0ftrt(2) = XYZE(3 )RCCEF( 1 ) - XYZE(1 )RCCEF(3 ) 
VNORM(3 ) = XYZE (1 )RCCEF(2) - XYZE (2 )RCCEF(1 ) 

RVEHMG s VXYZE (I) 2 + XYZE (2) 2 + XYZE (3) 2 

T3 = VNORM<R )R VEHMG 

T4 = VNORM ( 1 )XYZE(2 ) - VNORM( 2 )XYZE ( 1 ) 

BAR VEH = ATAN2(T3,T^) 

RCHMAG = VrCCEF(I) 2 + RCCEFC2) 2 + RCCEF(3) 2 
SINB = RT URN/RCHMAG 


Figure C-4.- EGRT-BV, bearing of vehicle. 
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EGRT-BVCHAC 


T5 = VN0RM(3 )RCHMAG 

T6 = RCCEF(2 )VN0RM(1 ) - VN0RMC2 )RCCEF ( 1 ) 
8ARCC = ATAN2(T5,T6) 


Figure C-5.- EGRT-EVCHAC, bearing to center of alinement circle. 
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EGRT-COSTHETA 



\ 

DO FOR I = 1,3 

/ 


> 


T7 = T7 + XYZE (I )RCCEF(I) 


CTHVC = T7/(RVEHMG RCHMAG) 


Figure C - 6 . - EGRT-COSTHETA, great circle arc. 
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Ht'FM.’ 1 

M_T-DW^F1__ 


GTHVB 'Vi. - CTHV<’- 

iTvwri = ivmvt: + u. ■ ■ th vc siNfv’ 

SBAROR = 1 NB/STli VV 

TEMP = CTVWP1 SBARCR 

TKMP = AMIN 1(1. ,AMAXKTEMP,-1.n 

A, 1 s AOOS(TEMP) 

TKMP = CTVWP1 

TEMP = AM IN 1 ( 1 . , AMA X 1 ( TKM P, - 1 , ) ) 
BVKWP1 = ACOS(TKMP)RCHMAG 
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EBftL 

_J 

-D.VNEE- 

! TEMP = SBARCR ! 
| TEMP = AMIN1 ( 1 . , AMAX1 , -1)) i 
! T8 = ASIN(TEMP) I 
! BARWP1 = BARVEH - SIGN(T8, RC(2)) 1 
! P = SIGN(1 . , RC ( 2 ) ) ! 
i A3 =0.5PI-A2+P# (RAZ - BARCC) | 

l 1 




\ 

IF A3 < - 0.003 > 

/ 

t t 

1 A3 s A3 ^ 2PI | 

i i 



l 

1 

DARC = A3 RTURN ! 

TRANG = CNMFS (DVEWP1 DARC - XNEP) i 

i 

i 




\ ! 

IF VE > VMIDPT \ THEN ! 

AND IPASS1 = 0 / | 

/ ! 

TRANGP = TRANG j 


Figure C-8.- EGRT-DVNEP, range- to- threshold point. 
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EGRT-DVNEP 


i 

i 

i 

i 

! 


1 

! [VE] r [TTD2IN]*- 1 [XYZED] 

1 

1 

I 

i 



FSI = AT AN 2 (VE (2) , VE(1)) 
DELAZ = PSI - BARWP1 


\ THEN 

IF DELAZ i > it V 
/ 





\ ! 



IF VE > VMIDPT 

\ THEN | DELAZ P = DELAZ 



AND IPASS1 = 0 

/ 




/ ! 



Figure C-9.- EGRT-DELAZ, azimuth error. 
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APPENDIX D 


IBM AUTOPILOT FLOW CHARTS 
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APPENDIX D 

IBM AUTOPILOT FLOW CHARTS 

The flow charts presented in this appendix represent IBM's response to the 
requirements set forth in this document. The autopilot model for the MCC 
will be programed by IBM from the flow charts contained in this appendix. 
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Figure D-l.- Continued. 
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Figure 11-1 Continued. 
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SUBTRACT OFF DEGREES TRAVELED TO GET 
TO MAXIMUM ROLL RATE (RRM) 

(RRM - RR) Cl IS TIME REQUIRED TO 
GO FRGM RR TO RRM. T2 IS TOTAL TIME 
REQUIRED TO REACH MAX. ROLL RATE (RRM) 
SECOND TERM IS TIME REQUIRED TO ROLL AT 
MAX. RATE TO SLOW DOWN POINT. RRM * Cl 
IS DE-ACCELERATION TIME. 


DIFF1 = 

DIFF ■ 

■ (C4 - C5) 

RR1 = 

RRM 


T2 

T1 + 1 

[RRM - RR) * Cl 

T3 

T2 + 1 

(DIFF! - C4 ) / RRM 

T4 

T3 + 1 

}RM * Cl 


IF DTIM > T2 


DOSEGMENT 

TIMCHK 


COMPUTE INCREASE IN ROLL RATE SUBTRACT 
OFF DEGREES TRAVELED WHILE INCREASING THE 
ROLL RATE. 

ELSE DTIM = DTIM - T1 I 

— RR = RR + RA * DTIM 

DIFF = DIFF + C5 - (RR * RR) * C2 


/ 


RETURN 

V 

J 


DOSEGMENT 

FINISH 


SEGMENT NORM 


Figure D-l.- Continued. 







THIS SEGMENT IS USED FOR TWO PURPOSES. 

1 . MOST COMMON USE IS WHEN DIFF IS SHORTER THAN 
DISTANCE TO SPEED ROLLING TO MAX (RRM) PLUS 
SLOWDOWN DEGREES (C4). HERE RR1 IS MAX VELOCIT 
REACHED AT START OF SLOWDOWN. 

2. A SECOND USE IS WHEN INITIAL ROLL RATE IS AWAY 
FROM ROLLC. HERE RR1 IS MAX ROLL RATE 
REACHED COMING BACK TO ROLLC AFTER STOPPING. 

T2 IS THE TIME AT RRl 

T4 IS TOTAL TIME TO REDUCE DIFF TO ZERO. 


ENTER 


SQRT (RA * (DIFF + C5)) 
T1 + (RR1 - RR) * Cl 
T2 

T2 + RR1 * Cl 


IF DTIM > T2 


DOSEGMENT 

TIMCHK 


ELSE DTIM = DTIM - T1 

RR = RR + RA * DTIM 
I DIFF = DIFF + C5 - (RR * RR) * C2 


DOSEGMENT 

FINISH 


SEGMENT UPDOWN 


Figure D-l .- Continued. 
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IF DTIM 
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Figure D-l.- Continued. 










SEGMENT FINISH 


Figure D-l Concluded. 
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